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Abstract. Conformal field theory has turned out to be a powerful tool to derive two- 
dimensional lattice models displaying fractional quantum Hall physics. So far most of 
the work has been for lattices with open boundary conditions in at least one of the two 
directions, but it is desirable to also be able to handle the case of periodic boundary 
conditions. Here, we take steps in this direction by deriving analytical expressions for a 
family of conformal field theory states on the torus that is closely related to the family 
of bosonic and fermionic Laughlin states. We compute how the states transform when 
a particle is moved around the torus and when the states are translated or rotated, and 
we provide numerical evidence in particular cases that the states become orthonormal 
up to a common factor for large lattices. We use these results to find the S'-matrix of the 
states, which turns out to be the same as for the continuum Laughlin states. Finally, 
we show that when the states are defined on a square lattice with suitable lattice 
spacing they practically coincide with the Laughlin states restricted to a lattice. 


Keywords: fractional QHE (theory), conformal field theory (theory), spin chains, 
ladders and planes (theory). 

1. Introduction 

In certain two-dimensional systems consisting of many interacting bosons or fermions, 
it is possible to have quasi-particles with unusual properties like fractional charge and 
an exchange statistics that is neither bosonic nor fermionic. The properties of such 
topological systems are very different from other types of matter, and they are therefore 
attracting a lot of attention. 

The quantum states appearing in connection with the fractional quantum Hall 
(FQH) effect constitute important examples of topological states, because they can 
be realized experimentally in semiconductor devices under demanding conditions, and 
because many of the states are believed to be described by simple analytical wave 
functions ^m\. The interest in topological systems has triggered the question whether 
FQH physics can occur in other settings, and in the last years there have been several 
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proposals for how to obtain FQH behaviour in lattice systems. This research allows 
for the investigation of new aspects of FQH physics and paves the way towards new 
possibilities for realizing the effect experimentally [IIH2n]. 

There are mainly two strategies to obtain lattice models displaying FQH physics. 
One of them is to engineer a band structure that is reminiscent of a Landau level, ensure 
that the band is partially hlled, and add interactions between the particles |2T1 - l2l] . The 
other strategy is to start from an analytical FQH wave function in the continuum, 
modify it to a lattice wave function, and then analyze its properties and use analytical 
tools to derive a Hamiltonian for which the state is the ground state [T^[251 - 136] . Both 
approaches have been shown numerically to be successful. 

It has long been known that there is a connection between FQH physics and 
conformal held theory (CFT) and that a number of continuum FQH states can be 
expressed as conformal blocks of certain CFT correlators jTj. This observation is 
also helpful for hnding analytical expressions of FQH states with periodic boundary 
conditions |9]. More recently, the CFT formulation has turned out to be very fruitful 
for constructing FQH models in lattice systems with analytical wave functions and 
corresponding parent Hamiltonians [T ^ l28 ] l2 ^lMl - l38] . 

A natural way to transform a wave function in the continuum into a lattice wave 
function is to keep the analytical expression for the state, but to let the state be a sum 
over all possible distributions of the particles on a lattice rather than an integral over 
all possible positions of the particles in the plane. Part of the trick used to derive the 
parent Hamiltonians in [28 ] l2 ^lMl - l35] l38]. however, is not just to do this, but to make 
an additional change. The Gaussian factor appearing in the continuum FQH wave 
functions is obtained by including a uniform background charge in the CFT correlator, 
and the additional change is to restrict the background charge to the same lattice as the 
particles. This has the further advantage that the CFT correlator takes a particularly 
simple form: for a lattice with N sites, it is the expectation value of a product on 
N CFT operators without any other factors. The analytical expression for the lattice 
states then differs slightly from the analytical expression for the continuum states, but 
at least as long as the distribution of the lattice sites is not too far from uniform, this 
does not change the physical properties of the states signihcantly. 

Apart from a study of the chiral spin liquid Kalmeyer-Laughlin state in [39], the 
investigations of lattice FQH behaviour obtained from CFT as described above has so 
far been concerned with systems having periodic boundary conditions in at most one 
direction, because the case of periodic boundary conditions in both directions is more 
complicated to handle in CFT. It is, however, very desirable to be able to investigate 
two-dimensional models with periodic boundary conditions in both directions, since the 
topology of the space on which a state is dehned plays an important role for topological 
systems. In particular, several computational methods have been developed in recent 
years to determine particular topological properties of quantum systems [lOflTF] . but 
many of these rely on the ability to dehne the investigated systems on surfaces with 
periodic, or more generally twisted, boundary conditions. Systems dehned on a torus 
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are also appealing because the lack of a boundary makes it easier to reliably investigate 
the bulk properties of the states. 

In the present paper, we take steps towards generalizing the previous work on lattice 
FQH physics obtained using CFT techniques to the torus. Our starting point is a CFT 
correlator that has previously [31] been used to derive lattice models with Laughlin-like 
ground states on the plane and cylinder. By evaluating the same correlator on the torus 
we obtain analytical expressions for a family of CFT states that is closely related to 
the family of bosonic and fermionic Laughlin states on the torus. We investigate the 
boundary conditions of these states when a single particle is moved around the torus, 
and we study how the states on particular lattices transform under collective rotation 
and translation operations. We provide numerical evidence in particular cases that the 
states become orthonormal up to a common factor for large lattices, and from this and 
the properties under rotation we conclude that the S'-matrix of the states is the one 
expected for the Laughlin states. The ij’th entry of the S'-matrix gives the phase factor 
acquired by the wave function when an anyon of type i is moved adiabatically in a closed 
path around an anyon of type j. The result hence shows that the braiding properties 
of the anyons are those of the Laughlin anyons. Finally, we show that, for the case of 
a square lattice, the analytical expressions of the states practically coincide with the 
analytical expressions of the corresponding Laughlin states on the torus. 


2. Lattice Laughlin states on the plane from CFT 


We hrst briefly recall the construction of lattice Laughlin states on the plane from CFT 
proposed in [31]. The starting point is the CFT correlator 

{(t)nAzi^zMn2{z2,Z2) ' ' ' , Z^))^, ( 1 ) 


where (...) denotes the vacuum expectation value, the subscript P stands for plane, and 


f’rij {Zj, Zj) 


. 9^7-1 

1 - 

e 




( 2 ) 


is a vertex operator [38]. Here, : ... : means normal ordering, g is a positive integer, 
Hj G {0,1}, Zj is a complex number, Zj is the complex conjugate of Zj, and ip{zj,Zj) is 
the held of a free, massless boson compactihed on a circle of radius R = y/g. 

CFT correlators can be decomposed [381139] into conformal blocks that are functions 
of the holomorphic coordinates Zj, but not of the antiholomorphic coordinates Zj. This 
decomposition allows us to write ([1]) in the form 

{(j)m{zi,Zi)(j)n2{z2,Z2) ■ ■ ■ (j)nr,{zN,ZN))p = cl f^ {Uj, Zj) ff {uj, Zj), (3) 

k 

where /I”(uj, Zj) are the conformal blocks and c^ are numbers that do not depend on Zj 
or Zj. 

On the plane, there is only one conformal block ff {nj, Zj), and we use this conformal 
block to dehne a wave function on a two-dimensional lattice as follows. We interpret 
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(Re(zj), Im(zj)) as the position in the plane of the j’th lattice site and Uj as the number 
of particles on the j’th lattice site. We then write 

I#) = #(ni,n2,... ,nAr)|ni,n2,.. ..um), (4) 

n\ ,n2v?^iv 

where 


V'p(ni, ^2, . . . , Uat) OC fiiUj, Zj) = 5nXn 

i<j 


^ ypiiTLj-ni-nj+llq 


The factor bn is dehned as 


5 


n 


1 for = iV/g, 

0 otherwise, 


(5) 

( 6 ) 


and it hxes the lattice hlling factor ^o 1/g. Since (jS]) does not £x the phase of 

the conformal block, we include an unspecihed phase factor Xn- Note that Xn cannot 
depend on the lattice positions due to the requirement that the conformal block is 
holomorphic, but it may depend on Uj. 

As shown in [3T], the factor ni<j(A — approaches the Gaussian factor 

of the Laughlin state with hlling factor 1/g up to a phase factor when iV —)■ oo if the 
lattice is regular (with the area per lattice site being the same for all lattice sites) and 
the boundary of the lattice is a circle. The phase factor can be transformed away if 
desired, and in practice, N ^ 100 is already enough to obtain very good agreement ED. 
We thus observe that ([5]) is practically the Laughlin state with hlling factor 1/g, except 
that the possible positions of the particles are restricted to a set of lattice sites. In |31] , 
it has also been shown that these states have the same topological entanglement entropy 
— ln(g)/2 as the continuous Laughlin states. 


3. Conformal blocks on the torus 

In this section, we derive the conformal blocks of the correlator ([T]) on the torus (the 
hnal result can be found in section l33|l . We first dehne the torus and the coordinates, we 
are using. Let uji and UJ 2 be two complex numbers. The modular parameter is dehned 
as r = and we assume that Im(r) > 0. The torus is then the parallelogram 

spanned by uJi and UJ 2 with periodic boundary conditions. In other words, if n and 
m are integers, and a and b are real numbers in the interval [0,1[, then we identify 
the point (n + a)ui + (m + b)u 2 with the point aui + buj 2 - Let 2 ; be a point within 
the parallelogram. We then dehne the scaled coordinate f = z/uji, which we shall use 
throughout the paper. 

We would like to determine 

(0ni(6,6)0n2(6,6) • • •0njv(^iV)^Ar))T = Ckfk{nj, fj) fk{nj, fj) , (7) 

k 

where T stands for torus. Our starting point for evaluating ([7]) is the following relation 
for the correlator of a product of vertex operators on the torus [50] (see also ED) 
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ww 


(p,p)gr 


( 8 ) 


On the left hand side of (IH]), holomorphic (antiholomorphic) part of 

the held = ^{^j) + massless boson and the numbers {uj, Uj) 

G r, where F is the lattice of the momenta consisting of the elements (p,p) with 


n mR 



n 



mR 


n,m ^ Z. 


(9) 


Here, R = y/q is the compactihcation radius mentioned above. For the correlator (j7|). 
the relevant values of {uj, Vj) are (—1/y^, —1/y/q) and {{q — 1)/ ■^/q, (g — 1)/y^), which 
are obtained for m = 0 and ?7, = —lor?7, = g — 1, respectively. 

On the right hand side of ([H]), = 1 if Ylj — 0 ^^^1 5y = t) otherwise, 


r7(r) = 

n=l 


is the Dedekind eta function. 






i<j 


( 11 ) 


and is the function obtained by taking the complex conjugate of Hq and 

replacing Ui with z/j and p with p. E is dehned as 


^(6-e„r) = 




where the Riemann theta function is given by 


ei{i.T)^d 


1/2 

1/2 




e 


a 

Vb\ 




i7rT(n+a)^+27ri(nH-a)(^H-fe) 


nGZ 


( 12 ) 


(13) 


The only part of ([H]) that does not immediately factorize into holomorphic and 
corresponding antiholomorphic parts is the factor 

exp I iTTrp^ + 27 rip fiUi - mfp^ - 27 ripy^ | 

(p,p)er \ i i ) 

= exp [ivrrg (n/g + m/2)^ + 27ri (n/g + m/2) ^/y^] 

n,m£'Z 

X exp [—infq {n/q — m/2)‘^ — 27ii {n/q — m/2) ^/^/q\ , (14) 


where we have dehned f = gX)i^*D- We shall consider the cases g even and g odd 
separately in the following. 


3.1. Even q 

When g is even, (n + mg/2) and {n — mg/2) are integers and differ by a multiple of g. 
Therefore, upon division by g, they give the same remainder 1. We can thus write 

p= ^{r + l/q), p= y/q{s + l/q), 
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r = 


n + mq/2 

S — 

n — mq/2 

(1 

5 

(1 


(15) 


Z = (n + mq/2, mod q) = {n — mg/2, mod g), 

where I G {0,1, ..g — 1} and r,s ^ 7^ without any restriction. That is, for a hxed r, the 
summation over s is over the whole set Z. Therefore, (ITTll becomes 

" ^^i^qr{r+lY+27ri{r+l)^^^-inqf{s+lY-2ni{s+^)^^ 


«G{0,l..g-l} rGZ 


which is a sum of factorized terms. 

3.2. Odd q 

When g is odd, {n + mq/2) and {n — mq/2) can take half integer values. In this case we 
split the sum over m in flTT|) into a sum over even and a sum over odd m. When m is 
even, we again write 

p=^(r + //g), p=^{s + l/q), 


r = 


n + mq/2 

s — 

n — mq/2 

Q 

5 

Q 


(17) 


Z = (n + mq/2, mod q) = {n — mq/2, mod g), 

but here we note that r and s are not unrestricted. Particularly, since m is even, we 
must have 


p — p 


= m = r — s = even. 


(18) 


Therefore, for a hxed r, say even (odd), the summation over s is over even (odd) integers. 
We can hence rewrite the sum in (ITTll over even m into 


i7rgr(r+i)^-i7rgf(s+i)%27ri(r+i)^-27ri(s+i)W f 1 + 


^i7r(r—s) 


E 

r,seZ 

where the last factor is 0 when r and s are of opposite parity and 1 otherwise. 
When m is odd, we write instead 

p= ^/q{r + l/q +1/2), p = y/q{s + l/q + 1/2), 
r = 


(19) 


n + (m — l)g/2 

s — 

n — [m — l)q/2 





( 20 ) 


Z = (n + (m — l)g/2, mod q) = {n — (m — l)q/2, mod g). 
In this case, we have 

p — p 




= m = r — s = odd. 


( 21 ) 


and so the summation over s has opposite parity to that over r. This allows us to rewrite 
the sum in flTdll over odd m into 


ijrgr(r+i + i)'^-i7rgf(s+i + i)%27ri(r+t + i)^-27ri(s+t + i)X fl 


r,sgZ 

The expression in flTT)) then equals the sum of (IT^ and 


' vV 


.( 22 ) 
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3.3. Expressions for the conformal blocks 

Collecting the above results, we conclude that the conformal blocks take the form 


fl,A{'^ji ^j) 


^nXr 

Vir) 


e 


l/q + aA 
bA 


iiiqui - 1 )> I n ~ 


quiTij—ni—nj+l/q 


',(23) 


i<j 


where = 1 if ^ • n, = N/q and = 0 otherwise, Xn is an undetermined phase factor 
that may depend on nj, but not on fj, and we have split the index k into two indices I 
and A, where I G {0,1,..., g — 1}. For q even, A can take only the value I, and for q 
odd, A G {I, II, III, IV}. The numbers a a and bA are defined as 


(24) 


~Q,iv = 1/2. There are hence q conformal blocks for q even and 4g conformal blocks 
for q odd. 


czi = 0, 

bi = 0, 

(for q even or odd). 

an = 0, 

bii = 1/2, 

(for q odd). 

am = 1/2, 

bill = 0, 

(for q odd). 

aiv = 1/2, 

biY = 1/2, 

(for q odd). 

e factors qj 

in (I7j) are all 1, 

and for q odd, we have 


3 . 4 . Phase factor 

Let us hnally discuss the relation between the phase factor Xn for the conformal blocks 
on the torus and the phase factor Xn for the conformal block on the plane. Locally the 
torus looks like the plane, and we can hence recover the conformal block on the plane 
by putting all the N lattice sites close together. Note that this does not affect Xn, since 
Xn does not depend on fj. In this case. 




and 


9 


- l),g^ 




9 


(0,gr) 


(25) 


(26) 


The right hand side of (|26|) is zero for a = 6 = 1/2, but is nonzero for the other values 
of a and b occurring in the centre-of-mass factors. In this limit, we hence have only one 
conformal block, satisfying 

/fc(e„ n,) OC (X Hiz, - ,^yn.n,-n.-n,Allq_ ^27) 

i<j i<j 

Comparing this result to ([5]), we conclude that it is natural to choose Xn to be the same 
on the torus and on the plane. In the following, we shall take 

x„ = (-yEiU-i)"., 

For q = 2 this choice ensures that the state is an SU(2) singlet 


( 28 ) 
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4. Wave functions from conformal blocks 


We use the conformal blocks on the torus fl23l) to dehne the wave functions 
\fjl/q+aA,hA) = (29) 




'4^l/q+aA,bAi^j 1 ^nXnd 


l/q + aA 
bA 


N 


- l),gr 


griinj—ni—rij 


. i=l 


i<j 


where V'z/g+aA,feA(^j) %) is shorthand notation for V'z/g+aA,feA( 6 )... ,fN,ni,..., un) and 
we have left out some constant factors that only affect the norm of the states. Note 
that \'ipi/q+aA,bA) is normalized. 

Since 


n E{fi - fq, = (-l)EAj-Ln, JJ JJ - fj, r)-"% (30) 

i<j i<j ijtj 

we can alternatively express (l29|) as 

\'4^l/q+aA,bA) ~ ^ V'pq+aA,bA('Cl5 • • • )'Cn;‘^1, • • • ;‘^m) 


f’l/q+aAMi^ly • • • Si, . . . , Em) — 0 


X 


n^(: 


l/q + aA 
bA 

M N 


M N 

i=l i=l 




T 




i<j 


i=l j=l 


Here, M = ~ is total number of particles, at is the operator that 

creates a particle at position S, | 0 ) is the vacuum state, and the hrst sum is over all 
possible conhgurations of the particles on the lattice, i.e., all {Si, S 2 ,..., Em} for which 
Ej G {.^ 1 , .^25 • • • 5 ^n} for all j and Sj 7 ^ Ej for all i and j. In addition. 


E(3,t) 


1 for S = 0 
E{E, t) otherwise 


(32) 


Since E{—E,t) = —E(E,t), we observe that the states (IMil describe bosons for q even 
and fermions for q odd. Note also that if we choose two of the coordinates Sj and Ej 
to be the same in (El]), then the amplitude of the term is zero because ii^( 0 ,r) = 0 . 
Independent of q, there is hence at most one particle on each site, which means that 
the bosons are hardcore bosons. This result reflects the property of the Laughlin states 
that the amplitude for two particles being at the same point is zero. 

In (J7|), we chose to order the vertex operators so that the one evaluated at the point 
fj was the j’th vertex operator in the correlator, and in fl29ll . we chose the ordering of 
the Uj’s in the ket to follow the ordering of the vertex operators in the correlator. We 
could instead have chosen the ordering such that the j’th vertex operator was operator 
number p{j) in the correlator, where p is some bijective map from { 1 , 2 , ...,A^} to 
{1,2,... ,iV}. In (|29|) . this would change Xn into ^ would change the 
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ordering in the ket, and it would change the product over i < j into the product over 
all i and j for which p{i) < p{j). In fl3T|) there is, however, no longer any reference to 
this ordering, and \'ipi/q+aA,bA) hence invariant under reordering. This means that it 
does not make a difference how we choose to label the lattice sites. 


5. Boundary conditions and fluxes 

In this section, we investigate the boundary conditions of the state fl3T|) when the fc’th 
particle is moved all the way around the torus and back to its original position, i.e. 
when Sfc —)■ Sfc + 1 or Sfc —)■ Sfc + r. We also discuss how to generalize the states to fulhl 
twisted boundary conditions. 


5.1. —)■ + 1 


When Sfc -)■ Sfc + 1, we have Utilizing 


e 


Ibi 


{^ + n,T) = 9 


a 

b + n 


(^, r) = e 


27r'ian 


lb 


If, t) for n G Z, 


we conclude that the centre-of-mass factor transforms as 


l/q + aA 
bA 




qr 


= 9 


l/q + aA 
bA 


e 


27ri(Z+ija^) 


The Jastrow factor transforms like 


n 

Kj 




^31 


r 




3i>k) 


n*«fc) - Sfc, r)-? UjOk) Ei^k - Sj, rf 

where we have used the property ± 1, r) = —E{f, r), and 

M N / . .\ -1 

Ey:^.^ -\- Oik \ / t\N 


Em V E{Zi-f^,r) 

The overall transformation is hence 


= (-!)". 


(33) 


(34) 


(35) 


(36) 


'fl/q+aA,bAi.^3i • • • 1 i ‘^ 1 ; • • • ; ^k—li ^k + 1 ; • • • ) ^m) 

= fjl/q+aAfiAi^i^ • • • Si, . . . , (37) 

This shows that for q even, the state is periodic, whereas for q odd it can be periodic 
{aA = 1/2) or antiperiodic {a a = 0). The phase factor obtained when taking a particle 
around the torus in the ui direction measures the flux through the hole of the torus. 
We thus observe that increasing I by one, increases the flux through the hole by one. 
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In this case, q Ei ^ q + d'r. We shall need the properties 


10 


if + pr, r) = 


'a 

lb 

a + n 

b 

E{f ± r, r) = 


a + p 
b 


pe 




(^,^) = (C,^), neZ, 


r . 


Now 

^^l/q + ttA 


+ qr 


= e 


l/q + ttA 
bA 


and 

n 

Kj 


Ely^i riSi}^ 

E{Ei-Ej,T) 


ni«fc) -^k-r, tY HiOfc) EiEk - Ej + T,r 


i(>k) - 


and 

M N 

nn 

i=i j=i 


n,«fc) - Sfc, r)9 nj(>fc) - Sj, ty 

— gEi(<fc)(-i^'r+i’r+27ri(Si-Sfc))ggJ]j.(>j.)(-i7rr-i7r-27ri(Sfc-SP)g 
_ ^-iTrTg(N/g-l)^iwg(N/g-l)^2wig J2i '^i-2iTigEk^-2TTig{N/g-l)’Ek 
_ ^-iiTT{N-g)^iiT{N-g)^2iTig J2i Sig-27ri7VSfc 


X -1 

E( Ei + rbik - ^j, t) 

E{Ei-fj,T) 

N 

— TT('g-i7f'r-i7r-27ri(Sfc-^j)\-l _ ^mTN+inN+2niNB.k-2niEYiij 


n« 

1=1 


(38) 

(39) 

(40) 


(41) 


(42) 


(43) 


Combining the factors, we get 

Yi/g+aA,bAi^i^ ■ ■ ■ ^^n,'E^i, ■ ■ ■, Sfc-i, Ek + r, Ek+i,..., Em) 

= Yi/g+aA,bAi^i^ ■ ■ ■ ,^n,'E.i, ..., (44) 

The state is hence periodic for q even, and it is periodic (6^ = 1/2) or antiperiodic 
{bA = 0) for q odd. We also observe that the flux through the tube of the torus is 
independent of 1 . 


5.3. Twisted boundary conditions 

From the above derivations it is clear that the states can be modihed to have twisted 
boundary conditions with twist angles 0i and 02 by modifying the centre-of-mass factor 
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i i / 

without modifying other parts of the wave functions. In this case 

t/’(/)l/(27r(j')+Z/g+aA,—</<2/(27r)+feyl ('Cl) • • • ) i ^1; • • • ; —1; “1“ 1; ^fc+l; • • • ) ^m) 

= 'C’</>i/(27rq)+Z/g+aA -?i2/(27r)+6A('Cl) • • • ,Cv,Si, . . . , (-1)^^^“^“^^ 

and 


0i/(27rg) + l/q + aA 
-02/(27r) + hA 


(45) 


(46) 


0(/)l/(27rq)+Z/ij+aA,—</<2/(27r)+fe^ ('Cl; • • • i ^Ni ‘^1; • • • ) ‘^fc—1) 4“ '^i • • • ) ‘^m) 

= 'f<t>i/(2TTq)+l/q+aA-4>2/(2TT)+hA (Cl) • • • ) Ca^, , . . . , . (47) 

A flux unit can be added through the hole or the tube of the torus by increasing 0i or 
02 by 27r. 


6. Transformation properties of the states 

Our next aim is to determine how the states transform under rotations and translations. 
The transformations O that we shall consider below can all be expressed in terms of 
a rearrangement operation d~^{j), where d is a bijective map from {l,2,...,iV} to 
{1, 2,..., iV}. After the transformation, the number of particles on the site at position 
fj is nd-i(j). We can, however, equally well state that the number of particles at site 
fd{j) is Uj. The transformed state is hence given by 

C^l'fl/q+aAybA) ~ ^ ^ 4’l/q+aA,bAi^d{j),nj)\ni, . . . ,n]\f). (48) 

ni,...,n7v 

For the case of rotations, we shall assume that r = i and that the lattice is 
invariant under a 90° rotation around the origin, i.e. that the sets {i.^i, i.^2, • • •, iCv} 
and {Cl) C2, • • •, Cvf} are the same. For the case of translations, we shall assume 
that the lattice is an x Ly square lattice with lattice constant c. In this case 

oji = cLa,, UJ 2 = icLy, and r = iLy/L^. We shall also use the auxiliary indices 

Xj G (0,1... La; — 1} and yj G (0,1... L^ — 1} that are dehned such that the index 

j = Xj + i/jLa; + 1. In this notation, the position fj of the j’th lattice site is 

Ci = [^j ~ {Lx — ^)/‘^ + i{yj — {Ly — l)/2)]/Lx if we choose the origin such that Cj = 0- 


6.1. Rotation by 90° 


A rotation by 90° is described by the transformation fd{j) = Ci/'^- Note also that 
r = i = — 1/r. Let us first compute 


9 


C -q 


T T 


= X^exp 


^ (n + a)^ + 27ri(n + a) (— + h 


T 


E 

fcez ' 


exp 


mq 


2TTikx - {x + ay + 27ri(a; + a) ( —Vh 


T 


c 


r 


da; 



Lattice Laughlin states on the torus from conformal field theory 


12 


./-exp 


-2mak + \— + h + k 

q 


(49) 


Now write k = qn + m, where m G {0,1,... g — 1}. Equation (I49|) becomes 

9-1 


e 


'T 'T i < ^ ^ \i iq 


m=0 n€Z 


-27ria(gn + m) + ivrrg ( n H-] + 


m + b\ ivr,^^ 


q 


qr 


+27ri,^ ( n + 


b + m 


q 


9-1 


m=0 ndh 


mrq ( n + 


m + ivr.^ 


( T- i.rc 2 1 ^ r "t+fe 

^—hS-Triab 


-e 


9 

m=0 


^ ^ 6 + m' 


H-h 27ri n + 


9 

—ga 


gr 

(^,gr). 


— go) + 27ria6 


A similar computation gives 
E 

and hence 




r r 




(50) 


(61) 


/f f . N griinj-ni-nj 


i<j 


i<j 




quirij—ni—rij 


N^\N N iTrg / ^ ^ ^ T—T 


qUinj—ni—Tij 


(52) 


i<i 


where we have used the properties — 0 — 0 ; which follow from the 

assumption that the lattice is invariant under a 90° rotation. 

Combining the factors and utilizing r = i and fl33p . we conclude that 

1 9-1 


'<Pl/q+aA,bA{id{j),nj) = i 2 . + 2 


V ^ m=0 


—27ri/m 


^ (Ol ^i) 


In particular, utilizing (13^ and (EH]), we obtain 
2 

— + —-— 


m=0 


9-1 




m=0 


9-1 




(53) 

(54) 
(65) 
(56) 


m=0 


9-1 


i^i+1 i(edO-),n,-) = ^57^ 

g" 2'2 g 2 2 


m=0 























Lattice Laughlin states on the torus from conformal field theory 


13 


where it is assumed in the last three equations that q is odd and TZim = e 27riZm /(}^^ 
are the elements of the unitary matrix 



1 1 
X 0 - 2 vri/'? 

X e-27ri((?-l)/g 


1 

Q-2ni{q-l)/q 


g-27ri(g-l)2/g 


(58) 


We observe that states with periodic (antiperiodic) boundary conditions in both 
directions transform into states with periodic (antiperiodic) boundary conditions in 
both directions, while states with periodic boundary conditions in one direction and 
antiperiodic in the other transform into states with antiperiodic boundary conditions in 
the one direction and periodic in the other. 


6.2. Translation by one lattice site in the x direction 

Translation by one lattice site in the x direction is described by the rearrangement 
transformation 


d-\j) = 


for Xj > 0 


J + 1 

j - {L^ - 1) for Xj = 0 
The corresponding change in lattice positions is 




for Xj > 0 


0 - 

fj + {L^ -l)/L^ for Xj = 0 


This leads to 


fi — fj for Xi,Xj > 0 or Xi,Xj = 0 
^d(i) - 1 for Xi = 0, Xj > 0 

fi — fj — 1 for Xi > 0, Xj = 0 


and 


- 1 ) = - 1) - Ly + quL, 

i i 

where = X]{i|xi=o}Utilizing fl5^ and ± l,r) = we obtain 


9 


'^id(i){qni-l),qT\ = 9 “ 


, gr e 


2ma{—Ly+qnL) 


and 


qriinj—ni—nj 


i<j 


i<j 


(59) 


(60) 


(61) 


(62) 


(63) 


rij—ni—nj , 


rij—ni—nj , 


. X ^ )-nL{N-Ly)-Ly {N/q-riL ) 

.I'^-gnL-LyN/q 


(64) 


i<j 
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Collecting the factors, we get 

In particular, 

'^l/qfl{id(j),nj) = nj) Q-‘^^'^Ly/ql^_^YnL+LyN/q^ 

V’pg,l/2(^d(i), rij) = i)l/q,l/2{^j, Uj) e-^^^lLy/q(^_YqnL+LyN/q^ j'g^^ 

i^i/q+i/2fl{^dij), rij) = 'fi/q+i/2,o{^j, Tij) (- 1 )^^'^ ( 68 ) 

f^l/q+l/ 2 ,l/ 2 {^dij),nj) = ^l^l/q+l/2,l/2{^j, %) /q ^_YLy{N/q-l) _ ^gg^ 


Note that (—= 1 for g even. We thus observe that the states |'0z/g+aA,6A) 
are eigenstates of the translation operator, whenever we choose periodic boundary 
conditions in the x direction. 


6.3. Translation by one lattice site in the y direction 

The calculation for this case is similar to the previous one. Now we have 


and 


leading to 


d (j) = 


^d(j) = 


j + Ly. for yj > 0 
j - {Ly - 1)L„ for yj = 0 

fj-i/Ly, for yj > 0 
fj + i(Ly - 1)/L„ for yj = 0 


fi - fj for yi, yj > 0 or yi, yj = 0 
6(0 - 6(i) = { ^i-^j + r for yi = 0, yj > 0 
ii Vi > 0 ’ Vj = 0 


and 


'^id{i)qni = + qrinB - Ly,/q), 

i i 

where ub = X]{j|yj=o}%- Using the property 


e 


{f + pr, t) = 


a + p 
b 


(6 r) for p G 


with p = Ub — Ly./q, we get 
l/q + Oa 


e 


^id(i){.qni - l),gr =9 


l/q + Oa + Ub - Ly./q 
bA 


(70) 


(71) 


(72) 


(73) 


(74) 


'^^i{qni - l),gr 
i / 

(75) 


y^^-i-irqT(nB-Lj:/q) -2wi{nB-Ly:/q)(J2i ^i{qni-l)+bA) 

Also, from fl40|) . we get 

6 JjE(^rf(i) 


i<j 


X 




quinj—ni—rij 


i<j 


(76) 
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After some computations, this gives the transformation 

In particular, 

f^l/q+aA^dij). rij) = ^P(l-L,)/g+aA,0{^j, rij) (-iy^B+NL,/q^ 

f’l/q+aA,l/2{U{j),nj) = f’il-L,)/q+aA,l/2{^j, ^j) (79) 


where q odd has been assumed in the last equation. The states with periodic boundary- 
conditions in the y direction and with L^/q integer are hence eigenstates of the 
transformation. 


7. Orthonormality and the S'-matrix 


In particular cases, it is possible to utilize the transformation properties derived in the 
previous section to analytically hnd linear combinations of the states that are orthogonal 
to each other, but it is not guaranteed in general that the states obtained from CFT 
are orthogonal. In £gure[T], we study the overlap between the states |'0z/g+aA,b7i) on an 
Lx L square lattice with periodic boundary conditions numerically for q = 2 and q = 3. 
We do this by rewriting the overlaps into an expectation value 

, , I , \ _ - nN)Ff^{ni, 

\Vl/q+aA,bA\xk/q+aA,bA) ~ / , 

where 


,nN) 


E 


pcm 
P I pA 


(ni,... ,n7v)|^ 


(80) 


Ffff{ni,...,nN) = 0 


l/q + Oa 
bA 


N 


- l),gr 


. i=l 


is the centre-of-mass factor and 


PA(ni,.. .,njy) = (5„^ |F'”(ni,... ,n«)|' 




quirij—ni—rij 


i<j 


(81) 


(82) 


can be interpreted as a probability distribution for the conhgurations rii,..., uat. This 
allows ns to compute {'f>i/q+aA,bA\'^k/q+aAfiA) to a constant that does not depend on 
I and k using a Metropolis Monte Carlo algorithm. The hgure shows that the states 
become orthogonal for large lattices. In addition, all the states have the same norm, so 
that they can be made orthonormal by multiplying by a common factor. Since the states 
\'f’i/q+aA,bA) have a well-defined flux through the hole and the tube of the torus, it follows 
that the S'-matrix for large lattices can be obtained directly from the transformation 
properties of the states under a 90° rotation |13]. This gives ns that the S'-matrix is 
the matrix 77 in fl58p (up to unimportant permutations of the rows and columns caused 
by the factors and This means that the phase factor acquired 

by the wave function when an anyon of type i is moved adiabatically in a closed path 
around an anyon of type j is exp(—27ripiPj/g), where pfc £ {0,1,..., g — 1} and Pk/q is 
the charge of an anyon of type k. We thus observe that the statistics of the anyons are 
as expected for the Laughlin states with hlling factor 1/g. 
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a) 1-6 


1.4 - 

1.2 - 


0. 

s 


0.6 


0.4 - 
0.2 - 
0 - 
- 0.2 


0 

e e 


q = 2 


b) 1-6r 


0.8 


0.2 - 

0 - 


g = 3 


0 (bilV’i) 


■—■ 0.6 

e 

0 (bilV’i) 

0 Re(('i/'i|b/fc>)L 7^ 


8 


0 Re(()/i;|)/'i)), 1 7^ k 

0 Im((V;ih/>*:)), 1 7^ k 




0 Im(('i/'j|V’i>)L 7^ fc 


0 


Figure 1. Overlap of the states on an L x L lattice for (a) q = 2 and (b) q = 3 

computed with Monte Carlo simulations, where lipi) = \'fi/q+aA,bA) with aj^ = Ba — 0 
for q even and qa — bA = ll2 for q odd. The overlaps in the figure are scaled such 
that ~ 9- states \tpi) become orthogonal and have the 

same norm for large L. They can hence be made orthonormal by multiplying all the q 
states with the same factor. A similar behaviour is found for the other possible choices 
of UA and Ba for q = 3. 


8. Connection to the Langhlin states on the torns 


We shall now compare the lattice states obtained from CFT to the continuum Laughlin 
states on a torus [ll|5l[52] in the Landau gauge (here we use the expression in [52] ) 


IV’^lt) OC 


d'^Zi ■ ■ ■ 


d^Z 


M 


l/q + {M-l)/2 
-{M-l)q/2 




g 


M 

X H e-(V2)(im(xpF \z,,...,Zm), le{0,l,...,q-l}. (83) 

i=l 

In this section, we shall often use the unsealed coordinates Zi and Zi instead of the scaled 
coordinates fi = Zi/ui and Sj = Zj/cai. As in previous sections, lower case letters refer 
to the positions of the lattice sites, while upper case letters refer only to the positions of 
the occupied lattice sites. In the Landau gauge, the coefficients of the wave functions in 
flH^ do not change when Z^ ^ Z^ + wi, and it changes by a Zj-dependent phase factor 
when Zi ^ Zi + W 2 (see e.g. [8] for details). 

Let us consider the case of an x Ly square lattice with — 0 lattice 

constant and choose uji = and U 2 = iy/^ Ly. In section 4.5 of [32], it has 

been shown that for this lattice and i G {1, 2,..., Af}, we have 

N 


IIAe 


Xi+yiLx+l 




-1 


OC 


-1) 


XiLy+Vi -■iT{yi-{Ly-l)/2)'^ 


(84) 


i=i 

(The proof in [39] is for N even, but it can easily be generalized to apply for both even 
and odd N.) Note in particular that f^i+yii^+i £ {'C 15 ^ 2 , ■ ■ ■, which means that the 
left and right hand side of fIM]) do not necessarily transform in the same way when the 
ith particle is moved by wi or W 2 . Under the assumption that Zi G {zi, Z 2 , ■ ■ ■, zn} for 
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all i, we can use 


to rewrite m into 

M 


l/q + aA 
bA 


M M 

X J], 

2=1 2=1 


E 


qZi 




qr 


n 4 


i<j 


Zi — Zj 

V V^La; ’ 


XiLy+Yi TTo-{l/2)(Im(Zi))2 




(85) 


where Xj and are defined such that Zi = \/^[Xi — {L^ — l)/2 + i{Yi — {Ly — l)/2)]. 
Utilizing (l2l)) and the properties fl3^ and fl3^ of the theta functions, we observe that 
the centre-of-mass factors in (l83|) and fl85|) are the same if 

{ I for q even 

I for g & X odd (86) 

IV for g odd & X even 

and Z —)■ / + g/2 mod g if both g and M are even. Note that X can only be odd for g 
odd, since the number of particles M = X/g must be an integer. The single particle 
phase factors be removed with a simple 

transformation if desired. The states in (155]) with the choice (18 6 p are hence practically 
the Laughlin states, except that the possible positions of the particles are restricted to 
a lattice. The states with A equal to II, III, or IV, or to I, II, and III for g odd have 
boundary conditions that differ by a sign in one or both of the two directions as can be 
seen from (1311) and (llTl) . 


9. Conclusion 

In summary, we have started from a CFT correlator of a product of vertex operators 
and used it to derive analytical expressions for a family of lattice states on the torus 
that are closely related to the bosonic and fermionic Laughlin states at filling factor 
1/g. We have computed various properties of the states and also shown how the states 
are related to the Laughlin states. The same CFT correlator has previously been used 
to derive Laughlin-like models on the plane and the cylinder 131], and the present work 
provides the corresponding wave functions on the torus. Finding the states on the torus 
is interesting because it opens up possibilities to use a number of tools to study the 
topological properties of the states. Here, we have used some of these to find the S- 
matrix, which shows that the braiding properties of the anyons are the same as for the 
Laughlin states. These results and previous investigations of the states in the cylinder 
and plane geometry [3Tl[3ll[36] show that the presence of the lattice and the small 
differences in analytical expressions between the CFT states and the Laughlin states do 
not change the topological properties of the states. In addition, the states on the torus 
are helpful for numerical studies of, e.g., the bulk properties of the states. 
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